Periodic-Orbit Theory of Anderson Localization on Graphs 
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We present the first quantum system where Anderson localization is completely described within 
periodic-orbit theory. The model is a quantum graph analogous to an a-periodic Kronig-Penney 
model in one dimension. The exact expression for the probability to return of an initially localized 
state is computed in terms of classical trajectories. It saturates to a finite value due to localization, 
while the diagonal approximation decays diffusively. Our theory is based on the identification 
of families of isometric orbits. The coherent periodic-orbit sums within these families, and the 
summation over all families are performed analytically using advanced combinatorial methods. 
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Anderson localization is a genuine quantum phe- 
nomenon. So far, attempts to study this effect within 
a semiclassical (periodic-orbit) theory seemed to be 
doomed to fail from the outset: It is not clear whether the 
leading semiclassical approximation for the amplitude as- 
sociated with a single classical orbit is sufficiently accu- 
rate. Even more seriously, there is no method available to 
add coherently the contributions from the exponentially 
large number of contributing orbits. Here, we address the 
second problem and develop a method to perform the co- 
herent periodic-orbit (PO) sums in a standard model — a 
quantum graph analogous to the Kronig-Penney model 
in ID — for which the PO theory is exact. For a list of 
references on the long history of graph models see Q. 

For investigating Anderson localization we consider the 
quantum return probability (RP). It is defined as the 
mean probability that a wave packet initially localized at 
a site is at the same site after a given time. We show that 
the long-time RP approaches a positive constant, which 
proves that the spectrum has a point-like component with 
normalizable eigenstates. The asymptotic RP is the in- 
verse participation ratio, which is a standard measure of 
the degree of localization. The RP can also be seen as 
2-point form factor of the local spectrum |2) . As such, it 
belongs to the class of quantities which can be expressed 
as double sums over PO's of the underlying classical dy- 
namics p]. Because of the exponential proliferation of 
the PO's in chaotic systems, the resulting sums are hard 
to perform. Consequently, most semiclassical approaches 
to spectral two-point correlations were restricted to the 
diagonal approximation where the interference between 
different PO's is neglected [§||@-|]. While this method 
is very successful for short-time correlations, it fails to re- 
produce long-time effects such as Anderson localization 
which are due to quantum interferences. In M the univer- 
sal long-time behavior of the form factor was related to 
universal classical action correlations between PO's of a 
chaotic system. A deeper understanding of how quantum 
universality is encoded in classical correlations is highly 



desirable but still lacking, despite some recent progress 
[§,||. This context is our motivation for developing the 
first PO theory of ID Anderson localization, although 
the phenomenon as such is well understood [|To|— [T^] . 

Quantum graphs exhibit both classical and quantum 
universal properties, which qualify them as model sys- 
tems in quantum chaos pj. The PO theory in graphs is 
exact. Hence, graphs are well suited for the study of PO 
correlations and their expected universal features. Re- 
cently, we reproduced the complete form factor of the 
circular unitary ensemble of 2 x 2 random matrices using 
a PO expansion in a simple quantum graph |fr5f . The 
new combinatorial tools developed there will be used to 
compute the quantum RP, by extending a method due 
to Dyson [fl9f for summing over orbits in a ID topology. 
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FIG. 1. Top: Quantum chain graph with vertices (dots) 
and bonds (line segments). The arrow shows the ini- 
tial state used to compute the RP. Bottom: The family 
[mo,mi,m2] = [1,2,2] of different but isometric orbits re- 
turning after topological time n = 10 to the vertex 0. 

We consider a quantum graph with a ID chain topol- 
ogy and use the intuitive notation as indicated in Fig. 1. 
The bond lengths Lj — Xj + i — Xj are disordered, i. e. 
pairwise rationally independent. On a bond j the general 
solution of the ID Schrodinger equation for wave number 
k is ^j(x) = Oj.-i- exp(+i k[x — Xj]) + aj.- exp(— i k[x — 
Xj + {\). The matching of solutions across a vertex is 
achieved in terms of a unitary scattering matrix o~j [rj , tj) 
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The phase factors on the r.h.s. account for the free mo- 
tion on the bonds 

pendent of k and are parametrized as Tj 



The matrices Oj will be assumed inde- 



iT]j. 
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sin r)j . A similar model was used e. g. in the analysis of an 
optical experiment demonstrating Anderson localization 
in the transmission of light through a disordered stack 
of transparent mirrors |16| . The vertex-scattering matri- 
ces Oj could also be computed by assuming (^-potentials 
at x = Xj, as in the Kronig- Penney model |l7|| . In the 
following we shall consider two situations: (i) random 
cr, where the m are independent random variables dis- 
tributed such that the corresponding transmission and 
reflection probabilities Tj = \tj\ 2 , Rj — \rj\ 2 are uni- 
formly distributed in the interval [0, 1], and (ii) constant 
cr, where Tj = T, R } ■ = R V j. 

At fixed k we consider the Hilbert space of coefficient 
vectors a = a,j tV where j goes over the vertices, and 
v = ±. We introduce a unitary operator 

Uj',v';j,v(k) — e lkLj (6+„i pSji j+„ tj | »+i + 



(2) 



The map defined by U(k) describes the time evolution of 
a coefficient vector. It is the natural object for the inves- 
tigation of the graph El. k is an eigenvalue of the graph 
Hamiltonian iff 1 is in the spectrum of U (k), that is when 
all the matching conditions ([j]) can be satisfied simulta- 
neously sl=U (fc)a. In terms of U we can characterize the 
degree of localization by calculating the averaged quan- 
tum RP 



p(n) = (|([/"K +;0 , + l 2 ) 



(3) 



as a function of the discrete topological time n of a state 
which was initially prepared as dj jV = 6j_o5 Vt+ (arrow in 
Fig. 1, top). The average (. . .) in (]|) is over a large A; 
interval, and for random cr also over the ensemble of a 
defined above. Due to the structure of U, V{n) ^ only 
for even n = 2m. 

The classical analogue of the quantum graph M is a 
Markovian random walk with vertex reflection and trans- 
mission probabilities which are equal to the quantum me- 
chanical ones defined above. A classical trajectory is en- 
coded by the sequence of traversed vertices {j u } which 
must be consistent with the connectivity of the graph. In 
our case v v = — j v = ±1- Given the initial vertex, a 
trajectory can also be identified uniquely by the sequence 
{«„}• Because of the probabilistic nature of the dynamics 
all trajectories are unstable. In general, the representa- 
tion of a quantum evolution operator in terms of classical 
trajectories involves a semiclassical approximation. Here 



it is exact and amounts simply to expanding the matrix 
products in (||). The result 
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~^A\ exp(i kC\) 



(4) 



can be interpreted in terms of the classical trajecto- 
ries introduced above. A runs over all trajectories con- 
tributing to the RP (||), i. e. all sequences {v v } (y — 
0, . . . , n — 1) with v v = 0. The traversed vertices are 



j v = Ylv'<u Vv ' ■ ^ or g ra ph s ; the concepts of returning 
trajectories and PO's coincide, hence (Q) is a PO sum. 
The length of an orbit is L\ = X)"=o Lj v , such that kC\ 
is the corresponding action in units of h. The amplitude 
A\ is the product of the transmission and reflection am- 
plitudes accumulated along the orbit A\ — Y\.u=i A\ >v 
with A\. v — tj u if f„+i = v v and A\. v = Tj u otherwise. 
It was shown in [Q that \A\\ plays the role of the sta- 
bility amplitude of the orbit, while the phase of A\ is 
equivalent to the Maslov index. 
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FIG. 2. Exact quantum and classical return probability 
V{2m) as given by Eqs. (§), @. 

Expanding |.. ,| in (^) we obtain a double sum over 
PO's of the type studied in j|. For short time, the inter- 
ference terms pertaining to two different PO's which are 
not related by an exact symmetry can be neglected due to 
the averaging applied. For constant a with T = R = 1/2 
this diagonal approximation simply amounts to counting 
the number of classical PO's with period n — 2m since 



all weights \A\ 



2 n coincide. Any such PO is repre- 



sented by a code word {v u } containing m letters + and — , 
respectively. According to the initial condition vi = +. 
Hence, each PO A in (Q) corresponds to a way of select- 
ing from the remaining n — 1 time steps those m with 
negative velocity, and we have 
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(triangles in Fig. 2). It is well known that the diagonal 
approximation yields the classical RP H. And indeed, 
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for long time (|J) shows the expected classical diffusion. 
The decay of the classical RP to corresponds to the 
obvious fact that there is no localization in the classical 
analogue of our model. 

In the following we will show that a constructive inter- 
ference between PO's with different number of reflections 
but equal lengths leads to a finite saturation value of the 
exact RP and consequently to localization Q. When 
(||) is expanded into a double sum, only pairs of orbits 
with equal lengths survive the k average. Hence, all rel- 
evant interferences are confined to families of isometric 
orbits. Suppose a PO A returning to bond (0, +) after 
2m collisions with vertices traverses the bonds (J, +) and 
(J,—) TUj times, respectively (X^ m j = m )- The length 
of this PO is C\ = 2 J2j m jLj- Thus, for rationally inde- 
pendent bond lengths Lj a family of isometric orbits con- 
tains all orbits sharing the set C = [mo, m±i, . . .]. In con- 
trast to the diagonal approximation, these are not only 
symmetry-related orbits. A simple example is shown in 
Fig. 1 (bottom) , but with increasing orbit length families 
can contain more and more PO's. Taking into account 
orbit families Eq. (||) becomes 



V(n = 2m) = V] 
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The outer sum is over the set T n of families, while the 
inner one is a coherent sum over the orbits belonging to 
a given family. The phase and amplitude of each PO 
depend on its itinerary. The phase equals the parity of 
half the number of reflections. Had one assumed that 
these phases are randomly distributed within a family, 
(p|) would again reduce to the diagonal approximation 
p0| . In contrast, the exact result derived below from 
Eq. (jq) is shown in Fig. 2. For both, constant and ran- 
dom cr, the RP saturates to a finite value indicating local- 
ization. Hence, quantum localization is due to a delicate 
and systematic imbalance between positive and negative 
terms within families. This can be regarded as a clas- 
sical correlation (deviation from a random distribution) 
of phases of PO's. The saturation value, i. e. the inverse 
participation ratio is expected to be inversely propor- 
tional to the localization length and the classical diffusion 
constant, see e. g. |(| and refs. therein. Indeed we find for 
constant a and a diffusion constant D = T/(l — T) 3> 1 
the relation \\m m ^ sx V{2m) = (Dir)- 1 Q. 

In the sequel we outline the analytical evaluation of 
Eq. (m , deferring a detailed expose to a subsequent pub- 
lication p0[ . A sum with a similar structure was calcu- 
lated by Dyson jl9| to study a disordered chain of har- 
monic oscillators. This system is analogous to our model, 
but Dyson computed a different quantity — the density of 
states. It is essentially different from the RP because the 
latter is a two-point correlation function. Dyson com- 
puted the traces of powers of a Hermitian matrix, and 
did not have to keep track of the phases encountered 



in the computation of diagonal elements of powers of a 
unitary matrix. Consequently the combinatorial argu- 
ments needed to evaluate (^) — though similar in spirit to 
Dyson's approach — are quite different from . 

To start, we consider a one-sided graph with j > 0, and 
Ro = 1, and perform the coherent sum over all orbits A 
in a given family. To keep track of the precise number of 
reflections for each orbit, we choose to specify the orbit 
by the segments with positive velocity v v = + . In the 
schematic representation of Fig. 1 (bottom) these are the 
arrows which point to the right. When the left pointing 
arrows are deleted from the bonds 1 and 2 in Fig. 1 (bot- 
tom) one obtains the arrow structures displayed in Fig. 3. 
The vertical displacement of the arrows provides the in- 
formation about the time ordering which is necessary to 
reconstruct the complete orbit from the right pointing 
arrows. The lower an arrow, the later it appears in the 
trajectory. The number of arrows to the left (right) of 
the vertex j is rrij—i (mj). An orbit can be completely 
specified by prescribing, at all vertices, in which order 
the adjacent arrows are traversed. Note that this local 
time order can be chosen independently for each vertex. 
As a consequence, the factorization 
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V([m , . . . ,m 6 _i]) = V h (m b _i) Vj(r 
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) (7) 



results, b denotes the rightmost vertex along the orbit, 
i. e. mj — Vj > b. The local contribution Vj(mj, mj—i) 
from a vertex will be determined in the following. 

A B C 



FIG. 3. This figure corresponds to vertex 2 in Fig. 1, which 
is approached twice from the left and from the right. The right 
pointing arrows can be arranged in 3 distinct ways, which 
correspond to the 3 decompositions of the number 2 into 2 
non-negative integers 2=1+1=0+2=2+0 (left to right). 

Each local time order at vertex j corresponds to one 
of ( mj ,^ mj 1 l 1 1 _1 ) possibilities to distribute mj identical 
objects over rrij—i sites, or — in other words — to decom- 
pose mj into mj-i non-negative summands (Fig. 3). The 
number of transmissions to the right at vertex j is given 
by the number of non-zero terms in this decomposition. 
We can express Vj(mj,m,j—i) as a sum taking into ac- 
count the number v of such transmissions, together with 
the associated amplitude. We obtain 
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rrij-i 
v 



v-l 



,2v mj+m,j-i—2v 



(8) 



The standard definition of binomial coefficients ensures 
that the sum can be taken over all integers v. The power 
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of tj in the expansion (||) is 2v since the number of trans- 
missions to the left and to the right is the same. (™?~x) 1S 
the number of decompositions of rrij into v positive sum- 
mands and C 71 '^ 1 ) corresponds to selecting the v non- 
zero terms from all m,j—i summands. The sum in (Q) is a 
special Kravtchouk polynomial — a well-known object in 
combinatorics ]2l]fL5| ]. 

With (Q), (Ji) we can perform the summation over all 
orbits A of a given family in (||). However, there are 
Ylb=i (T-i 1 ) = 2 m_1 decompositions of an integer m 
into positive summands rrij, i. e. the number of fam- 
ilies grows exponentially with time. To sum over all 
families we derive a recursion relation which dramat- 
ically reduces the number of terms needed. Group- 
ing in (^) all terms according to the number of traver- 
sal of the first two bonds too and toi, i. e. defin- 
ing B m (m ,mi) = J2m 2 r P([m Q ,m 1 ,m 2 , ...}), we find 
B m (s + l,t) = V(t,s + i)/V(t,s)B ro _i(s,t) for s > 
and B m (l, t) = V(t, 1) J2T=i B m-i(t, s). Summing over 
the second argument of B m (s, t) we obtain the combined 
contribution V m (s) from all orbits with period 2m which 
traverse the initial bond exactly s times before finally 
returning to it. The corresponding recursion relation is 



( pi] ) simplifies to 7^ av ^(mi,_i) = l/(mb_i + 1). Conse- 
quently we can initialize the recursion (|^) by Vm^ (m) = 
1/(to+1), which results in Vm V \s) = s/(m 2 + m). Sum- 
ming with respect to s we find that the RP of the one- 
sided graph is constant T' os (2m) = 1/2 for to > 0. For 
the unrestricted graph and m — > oo the double sum in 
( ]l0| ) can be approximated by a double integral yielding 
lirriTO^oo V(2m) — 7r 2 /3— 3 which is indeed the saturation 
value of the top curve in Fig. 2. 

In summary, we have shown that Anderson localiza- 
tion can be reproduced from PO theory only when clas- 
sical correlations are properly taken into account. Here, 
these correlations show up as exact isometries of families 
of PO's, and the non-random distribution of the phases 
within a family. The failure of the diagonal approxima- 
tion in a model where PO theory is exact, identifies the 
neglect of action correlations as the primary source of er- 
rors in semiclassical theories of localization and related 
problems. 

Support by the Minerva Center for Nonlinear Physics 
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V m (t) = v ( s > t)V m -t(s) (1 < t < m) 



(9) 
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The recursion is initialized using the elements V m (m), 
which are due to a single PO: the orbit bouncing to times 
between the vertices j = and j = 1. The RP for the 
one-sided graph is now V(2m) = J2T=i ^m(s)- 

For the unrestricted graph a returning trajectory can 
be composed of simple loops to the right and to the left 
from the initial bond. Both groups can be described by 
the results for the one-sided problem. Using similar ar- 
guments as in the derivation of (0) we find 



V(2m) = 



m—m r -\-s 



(10) 



We were able to solve the recursion (JlO) analytically for 
random a. After expanding | . . . | 2 in (|S) we get 

1 /' 7r / 2 

Vj aV \ m 3, m 3-i) = 5 / sin27?jd% Y,[iu,.ni, ,.//,;• 



rrij , , (— \y+ v ' frrij-i + rrij 

i + rrij + 1 z — ^ vv 1 \ v + v' 



rrij-i — 1\ / rrij — 1\ / mj-i — 1\ ( rrij — 1 



v - 1 



v- 1 



v 1 - 1 



v 1 - 1 
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. . (11) 

(rnj_i + rrij — 1)(to./_i + m,j)[m,j-\ + rrij + 1) 

The last equality follows from an identity which we 
proved previously p5| using recent developments in com- 
binatorial theory J22|. For the outmost vertex on an orbit 
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